Direct Experimental Evidence for the Hybridization of Organic Molecular Orbitals 
with Substrate States at Interfaces: PTCDA on Silver 
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We demonstrate the application of orbital fc-space tomography for the analysis of the bonding 
occurring at metal-organic interfaces. Using angle-resolved photoelectron spectroscopy (ARPES), 
we probe the spatial structure of the highest occupied molecular orbital (HOMO) and the former 
lowest unoccupied molecular orbital (LUMO) of one monolayer 3,4,9, 10-perylene-tetracarboxylic- 
dianhydride (PTCDA) on Ag(110) and (111) surfaces and in particular the influence of the hy- 
bridization between the orbitals and the electronic states of the substrate. We are able to quantify 
and localize the substrate contribution to the LUMO and thus prove the metal-molecule hybrid 
character of this complex state. 



The adsorption of 7r-conjugated molecules on metals 
has attracted considerable interest in surface science due 
to its importance for the understanding of contacts in 
organic electronic devices. PTCDA (3,4,9,10-perylene- 
tetracarboxylic-dianhydride) serves as an archetypical 
molecule in this field [IH3- In particular, its behavior 
in molecular monolayers on Ag surfaces has been studied 
intensely as a model system for chemisorptive molecule- 
metal interaction Rich experimental data exist on 
the geometric structure of these interfaces [BHU , but the 
bonding distances differ quantitatively from theoretical 
predictions 0, Moreover, spectroscopic techniques 
show a charge transfer from the metal into the lowest un- 
occupied molecular orbital of the free molecule (LUMO), 
but the underlying mechanism remains debated. Liter- 
ature mentions the involvement of several molecular or- 
bitals Q, Ag s- and d- bands [ll[ and additional local 
bonds forming between the carboxylic O atoms and the 
Ag surface [1] . Furthermore the Shockley surface state is 
affected by the PTCDA adsorption and may play a role 
in the bonding as well 12, 13j]. A thorough understand- 
ing of the interaction at the interface must be based on 
the correct description of the valence orbitals. However, 
information about the spatial distribution of molecular 
orbitals (MOs) is not easily obtainable. So far scanning 
tunneling microscopy (STM) is the most direct approach 
to image surface charge distributions with sub-molecular 
resolution 0, Q . Another experimental method is angle- 
resolved photoelectron spectroscopy (ARPES) [3, llq . 
since the angular intensity distribution pattern of or- 
ganic molecules is linked to the wave function geome- 
try of the respective MOs 1|| 17 1. However, interpreta- 



tion of ARPES data including final-state scattering [18 1 
based on atomic orbitals is challenging [la. A new ap- 
proach has been described very recently [19J which di- 
rectly links photoemission intensity and initial state wave 
function, enabling energy resolved tomographic imaging 
of the electron distribution in fc-space. This allows a di- 
rect comparison of experimental and theoretical charge 



distribution and a probing of the spatial structure of the 
bonding orbitals. 

In this Letter, we present photoemission data of the 
two highest occupied interface states of PTCDA on 
Ag(110) and Ag(lll) surfaces and compare our data to 
free-molecule orbitals as calculated by density functional 
theory (DFT). We demonstrate the imaging of the re- 
spective orbital structure and its correspondence to the 
free- molecule HOMO/LUMO. Note here that we retain 
these terms also for the respective orbitals of the ad- 
sorbed molecule for simplicity, although the LUMO is 
not unoccupied anymore (while the HOMO becomes the 
second occupied MO). Good quantitative agreement is 
achieved for the spatial distribution of the (non-bonding) 
HOMO and DFT. The now occupied LUMO deviates 
from its free-molecule form due to substrate induced 
modifications. They take the form of an s-like contri- 
bution to the orbital localized laterally at the center of 
the molecule. Our interpretation of the electronic in- 
teraction is compared to quantum-chemical calculations 
of PTCDA adsorbed on Ag(110) and DFT results 
for PTCDA on Ag(lll) [9(. The experiments were per- 
formed in a UHV setup composed of an organic molec- 
ular beam epitaxy chamber for sample preparation and 
a spectrometer chamber for the ARPES measurements 
(base pressure below 10 -10 mbar). The Ag substrates 
have been prepared by standard sputtering and anneal- 
ing cycles as described elsewhere [2l|. PTCDA, puri- 
fied by triple-sublimation, was evaporated from a Knud- 
sen cell at rates of approximately 0.2 ML per minute 
onto the substrate which was kept at room temperature. 
The monolayer coverage was controlled by monitoring 
the continuous quenching of the Shockley type surface 
state [lH and the characteristic PTCDA valence band 
spectra of the chemisorbed monolayer jij. All ARPES 
measurements were performed at sample temperatures 
T R3 80 K with a high-resolution photoelectron analyzer 
(Scienta R4000) in combination with a monochromatized 
VUV lamp using He I Q radiation (hv — 21.23 eV) at an 
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FIG. 1: Right column: Experimental fc^-dependent ARPES 
intensity of the HOMO (top) and LUMO (bottom) of PTCDA 
on Ag(110) recorded with HeI Q . Left column: Correspond- 
ing DFT calculations of the free molecule at the respective 
photoelectron energies. 



energy resolution of AE = 5 meV. The angle resolved 
mode of the analyzer covers a parallel detection range of 
up to ± 15° with a resolution of « 0.3°. An additional 
rotation of the sample was used to allow 2D fc-space map- 
ping. Molecular orbitals of an isolated PTCDA molecule 
were calculated within the framework of DFT using the 
ABINIT software package [23| . Norm-conserving pseudo- 
potentials with a cut-off of 50 Ryd and a generalized gra- 
dient approximation for the exchange-correlation energy 
and potential have been used. The calculated 2D- ARPES 
intensity maps are obtained from the Fourier transform 
of the HOMO and LUMO orbitals at the appropriate 
kinetic electron energies. Details of this approach are 
described elsewhere [19(. 

The monolayer of PTCDA on Ag(110) is a single- 
domain structure with one molecule per unit cell. The 
identity of all molecules in terms of adsorption site and 
absolute orientation allows the effective mapping of MOs 
without the superposition of additional orientations. The 
molecule grows on Ag(110) oriented along the [001]-axis 
of the substrate 24( . We will refer to this direction in mo- 
mentum space as k x , to [110] as k y . Fig. 1 compares the 
experimental ARPES intensity (right column) recorded 
at binding energies of 0.8 eV and 1.9 eV with respect to 
the Fermi level to the corresponding calculations for the 
valence MOs (left column). The shown data set was com- 
pleted from measurements covering one quadrant using 
the symmetry axis of the system |24| . From the overall 




FIG. 2: Top: ARPES intensity vs. k y at k x = for 1 ML 
PTCDA/Ag(110), EDC for indicated k y on the right. Bot- 
tom: MDC deduced from the MOs energy windows (red/blue 
dotted line, dashed line marks respective zero) and at inter- 
mediate energies (gray curves). The black lines indicate the 
dispersing substrate bands. The solid red curve shows the 
theoretical LUMO intensity (offset from zero for better com- 
parison) . 



similarity we conclude that the measured orbitals corre- 
spond to the HOMO and the LUMO respectively. How- 
ever, the experimental data exhibits additional photoc- 
mission intensity features which do not appear in theory. 
They are mainly caused by the sp-bands of the Ag sub- 
strate. Fig. 2 (top) shows an energy-resolved ARPES 
spectrum along the fc^-direction at k x = through the 
momentum map in Fig. 1. The data allows a clear dis- 
tinction between photoemission intensity from the MOs 
and the strongly dispersing substrate bands. Selected 
energy distribution curves (EDC) are displayed on the 
right. Fig. 2 (bottom) depicts the energy-integrated k y - 
depcndcnt photoemission intensity of the LUMO (red) 
and the HOMO (blue). Additional line scans (momen- 
tum distribution curves (MDC)) at intermediate binding 
energies (gray lines) are added to visualize the contribu- 
tion of the Ag bands, black lines guide the eye along the 
dispersion. The momentum distribution of the HOMO 
signal is basically featureless, in agreement with theory 
(cf. Fig. 1). The LUMO signal also generally resembles 
the free-molecule calculations. However, there is an ad- 
ditional intensity maximum at normal emission (ku = 0). 
Note that this is neither present in the free molecule 
nor accidentally caused by substrate bands. In addi- 
tion, the position of the dominating emission maxima 
at k y w 1.6 A" 1 is shifted to smaller absolute fc^-values. 



3 




FIG. 3: Top: ARPES intensity vs k\\ for 1 ML 
PTCDA/Ag(lll), EDC for indicated k\\ on the right. Bot- 
tom: MDC of the LUMO (red curve) and HOMO (blue 
curve) intensity, respective energy windows marked above, 
zero indicated by dashed line. The solid lines correspond 
to the calculated DFT intensities (offset from zero). Gray 
shaded area models the normal emission intensity as Gaus- 
sian (FWHM = 0.42 A" 1 ) 



Interpreting our ARPES results for PTCDA on Ag(110), 
we find the HOMO generally unaltered by the adsorption, 
therefore we rule out significant substrate admixture to 
this orbital, the HOMO is apparently not involved in the 
bonding at the interface. The LUMO generally resembles 
the free-molecular calculations. However, the shift of the 
characteristic main intensity maxima to lower k y values 
indicates a distortion of the orbital structure in this direc- 
tion (i.e. the short molecular axis). Moreover, additional 
intensity appears around fen =0. In real space this would 
correspond to a laterally node-free charge distribution at 
the center of the molecule. We therefore suggest a bond- 
ing of the molecule due to LUMO hybridization, taking 
the form of an s-like substrate contribution to the MO 
at the PTCDA perylene core. A comparison with recent 
quantum-chemical calculations reconfirms the observed 
negligible interaction of the HOMO 0. For the LUMO, 
a distortion along the short molecular axis and additional 
substrate-localized lobes are visible in the calculated or- 
bital structure [2fjj, matching our experimental findings. 
However, the calculations show a node (and therefore a 
parity change) of the hybrid orbital along the k x axis, a 
prediction which is not compatible with the experimen- 
tally observed normal emission intensity. 

In contrast to the (110) surface, PTCDA on Ag(lll) 
forms a multi-domain herringbone structure. Six sym- 



metry equivalent domains with two molecules per unit 
cell [3j result in twelve different molecule orientations. 
Therefore, measuring the polar-angle dependence of the 
PE intensity effectively means probing the averaged fen- 
dependence of the MOs. While this reduces the 2D-maps 
into one-dimensional data sets it gives two advantages: 
On the one hand a more complex geometry of the or- 
ganic overlayer prevents an overlap with intense sp-bands 
due to backfolding 22]. On the other hand one mea- 
sures all data points at an identical geometry between 
photon source and the emitted electrons, allowing for 
a strict quantitative analysis of the MDC in terms of 
electron distribution of the initial MO, since the modu- 
lating polarization factor A • k remains constant [l9j ]. 
Fig. 3 presents energy and fen resolved photoemission in- 
tensity for a monolayer of PTCDA/Ag(lll) (top, data 
symmetrized to fen = 0). The MDCs of the background 
corrected LUMO and HOMO (red and blue energy win- 
dow) are compared to the free-molecule DFT results in 
the bottom panel. The theoretical data was averaged 
over the azimuthal angle (emulating the multiple do- 
mains in the experiment) and offset from zero to allow 
for a constant experimental intensity. The HOMO is 
again characterized by a single, dominating maximum. 
The deviation for higher fc-values is probably caused 
by the grazing incidence of the light and high electron 
emission angle (O « 70°), setting an upper k limit for 
quantitatively reliable measurements in our experimental 
setup. Otherwise, the agreement between free-molecule 
theory and experiment is excellent. From this we con- 
clude an undistorted HOMO geometry and a weakly in- 
teracting character of this MO on Ag(lll), analogous 
to the Ag(110) substrate. For the LUMO, a slight in- 
ward shift of the dominating maxima is observed, as is 
again ARPES intensity at fen — 0. From this we deduce 
a similar bonding mechanism as described for Ag(110) 
above. There is, however, also a distinctive new peak 
at fen s» ± 0.9 A -1 , neither present in the DFT calcu- 
lations nor a feature of the substrate. It corresponds 
to a real-space periodicity of 7 A in the LUMO struc- 
ture. This value is almost twice as large as the typical 
node-distance within the undistorted MO (causing the 
dominating ARPES intensity maxima) and comparable 
to the lateral PTCDA dimensions (9.2 x 14.2 A) [3]. This 
length scale hints at a charge redistribution over the sin- 
gle PTCDA molecule, possibly due to inter-molecular in- 
teraction. On Ag(lll) PTCDA arranges in a herringbone 
structure with a nearest-neighbor distance below the van 
der Waals radii of the molecules, caused by electrostatic 
forces acting between the negatively^charged anhydride 
group and the aromatic core region Q . Such an interac- 
tion could cause the observed modifications in a strongly 
delocalized bonding orbital like the LUMO. Another in- 
dication is the absence of the considered ARPES feature 
in the previously discussed data for the Ag(110) surface, 
on which the PTCDA molecules are distinctly separated 
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from each other on specific adsorption sites. In order to 
quantitatively analyze the additional normal-emission in- 
tensity, we approximate it by a Gaussian as displayed by 
the shaded area in Fig. 3. The real-space correspondence 
is a node-free charge spread laterally (FWHM w 7 A) 
over the central part of the PTCDA perylene core. A 
first estimate assigns « 10 % of the total PE intensity 
to the area around fc|| =0, a value which may pro- 
vide a measure for the substrate contribution to the 
PTCDA-Ag hybrid state. Quantum-chemical calcula- 
tions for PTCDA/Ag(110) state a comparable (16 %) 
substrate contribution [ioj ], A detailed comparison with 
theoretical calculations for PTCDA on Ag(lll) is not 
directly possible, there are only projected charge densi- 
ties for PTCDA/Ag(lll) obtained by DFT calculations 
0. They describe the lateral distribution between the 
molecule and the substrate as an undistorted LUMO with 
additional charge accumulated around the central carbon 
ring of the molecule This description resembles the 
electron probability density that would result from the 
orbital geometries deduced from our data. 

In conclusion, we demonstrate the potential of orbital 
tomography by ARPES for the comprehensive under- 
standing of the fundamentals of the bonding of large or- 
ganic molecules to surfaces. On the particular example 
of the PTCDA molecule, chemisorbed on Ag surfaces, 
we show that the HOMO orbital does not show a signif- 
icant substrate admixture. It is therefore not decisively 
involved in the molecule-substrate bonding. In contrast, 
the LUMO features a significant distortion of the orbital 
spatial structure compared to the free-molecule and a 
prominent admixture of substrate states. A first anal- 
ysis estimates the substrate contribution to the LUMO 
as sa 10%, localized in a node-free charge accumulation 
at the central carbon ring of the perylene core. Our re- 
sults are in good agreement with quantum-chemical cal- 
culations for the MOs of PTCDA on Ag(110) and DFT 
calculations for PTCDA on Ag(lll), with the exception 
of the nodal structure of the LUMO hybridization on 
Ag(110). 
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